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Flow rate of polygonal grains through a bottleneck: Interplay between 

shape and size 

Ezequiel GoldbergP C. Manuel CarlevaroP^I^ Luis A. PugnalonPEt 


We report two-dimensional simulations of circular and polygonal grains passing through an 
aperture at the bottom of a silo. The mass flow rate for regular polygons is lower than for 
disks as observed by other authors. We show that both the exit velocity of the grains and 
the packing fraction are lower for polygons, which leads to the reduced flow rate. We point 
out the importance of the criteria used to define when two objects of different shape are 
considered to be of the same size. Depending on this criteria, the mass flow rate may vary 
significantly for some polygons. Moreover, the particle flow rate is non-trivially related to 
a combination of mass flow rate, particle shape and particle size. For some polygons, the 
particle flow rate may be lower or higher than that of the corresponding disks depending 
on the size comparison criteria. 
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I. Introduction 

The flowing properties of granular materials are key 
in a myriad of industrial processes. In this respect, 
the discharge of a silo or a hopper constitutes the 
most ubiquitous operation. Considerable attention 
has been paid to silo discharge in the literature. 
Since the initial works in the ’60 [HE] a num¬ 
ber of studies have reproduced the essential aspects 
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highlighted in the pioneering works. These are re¬ 
lated to the little or no effect played by the height 
of the granular column (see however a recent ex¬ 
perimental report on submerged silos 0): and the 
5/2 power dependency of the flow rate with the 
size of the aperture (see however a study over a 
wide range of apertures 0)- The irrelevance of the 
column height for the flow rate has been recently 
shown not to be connected with the Janssen’s effect 
mm- Less attention has been put in the mate¬ 
rial properties and particle shape. In general, ma¬ 
terial properties seem to have a marginal effect on 
flow rate laiiiiiQ]. However, the flow pattern may 
change from funnel flow to mass flow to mixed flow 
as friction is varied either for the particle-particle 
interaction or the particle-wall interaction liiin]- 

Particle shape is also a factor that has been con¬ 
sidered by a number of authors. However, unlike 
material properties, shape is not defined by a fi¬ 
nite number of parameters. For a particular basic 
shape, for instance rectangles, one can assess the ef- 
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feet of the few parameters that describe that shape 
(e.g., long side and short side). However, this can¬ 
not tell if a different basic shape (e.g., ellipses) will 
yield different results. As a consequence, the effect 
of shape has to be surveyed necessarily in a limited 
manner and the selection of what basic shapes to 
consider has to be led by the application of inter¬ 
est or by practical issues such as the limitations in 
fabrication of some shapes or their computer mod¬ 
eling. 

Studies on circular grains and on supercircles 
(rounded squares) in 2D show a significantly lower 
mass flow rate for the supercircles HU. This au¬ 
thors, compare particles of different shape that 
have the same larger diameter. We will see that 
this selection of matching size may bias the conclu¬ 
sions. A number of ulterior works have shown the 
same overall trend that mass flow rate is reduced in 
comparison with rounded particles if they present 
sharp edges [Hllallll]. A contrasting result was 
found by some works in which elongated rounded 
particles flow faster than circles (2D) or spheres 
(3D); however, they considered frictionless grains 
in one occasion m and in the attempt to validate 
with experiments they found difficulties in produc¬ 
ing particles of same mass and different shape |16j . 
One particularly interesting exception to the re¬ 
duced flowability of sharp particles is the case of 
tetrahedrons that seem to flow faster than spheres 
[14j . Unfortunately, this study compared spheres 
and tetrahedrons of different mass/volume, and the 
conclusion cannot be warranted to be general. 

There is a significant point generally overlooked 
in many of these studies. Since the relative size 
between the particle and the aperture is the main 
factor affecting the flow rate, comparison with dif¬ 
ferent shapes has to be done by fixing the particle 
size. For example, some workers have claimed that 
hexahedrons and ellipsoids flow faster than spheres 
based on experiments with different seeds, without 
paying much attention to the fact that the seeds 
with non-spherical shape used were also smaller in 
size m- The problem was recognized in other 
works when comparing tetrahedrons with spheres 
M- There is an arbitrary decision to make as to 
what defines if two particles of different shape are 
of the same size. The canonical solution to this is 
to define that particles must have the same mass to 
be considered of the same size. Notice however that 
some studies use a characteristic length of the par¬ 


ticles (e.g., maximum feret diameter) [TTJ [131 HI] • 
Since the shape is being changed, keeping a charac¬ 
teristic length constant will lead to different masses. 
There is no correct selection in this respect. The 
parameter that defines the effective size has to be 
chosen according to the particular application, for 
example. We will show that this choice has a sig¬ 
nificant impact on the mass flow rate observed in 
particles of different shapes. 

An additional aspect to the flow of grains is 
the actual number of particles outpouring per unit 
time. Although for a given shape this is propor¬ 
tional to the mass flow rate, in comparing results 
for different shapes the particle flow rate has a 
different proportionality constant for each shape. 
Moreover, the criteria used to define if particles of 
same size are compared will strongly affect the par¬ 
ticle flow rate. We will present results that show 
that to deliver a fixed number of particles per unit 
time, the use of different particle shapes requires 
the increase of the aperture in some cases and the 
reduction in others. 

II. Simulation 

We simulate a 2D silo by means of a discrete ele¬ 
ment method (DEM) using the Box2D library [TB] . 
This package uses a constraint solver. At each time 
step of the dynamics, a series of iterations (typi¬ 
cally 20) are used to resolve penetrations between 
bodies through a Lagrange multiplier scheme m- 
After detecting overlaps, the inelastic collision at 
each contact is solved and new linear and angu¬ 
lar velocities are assigned. The equations of mo¬ 
tion are integrated through a symplectic Euler al¬ 
gorithm. The time step St used to integrate the 
equations of motion is 0.025 d/g; with d the diam¬ 
eter of our circular grains and g the acceleration of 
gravity. Solid friction is also handled by means of 
a Lagrange multiplier scheme that implements the 
Coulomb criterion. Box2D is highly efficient when 
handling complex bodies such as polygons. 

The silo consist of a box 20d wide and 200d 
high. Gravity acts in the negative vertical direc¬ 
tion. 2000 grains are typically used in each dis¬ 
charge, and these are not reinjected after they leave 
the silo. After grains initially placed at random 
without overlaps settle, an aperture opened at the 
center of the flat base of the silo allows the grains 
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Figure 1: Snapshots of the simulations for disks, 
triangles and squares. Only part of each system is 
shown. 

to discharge. For small apertures, clogging may be 
observed. In such cases the arrested discharge is 
discarded from the analysis. Only those discharges 
that did not present clogging are considered. Some 
material may be left over at the stagnant zones on 
each side of the orifice. 

Grains of different shapes corresponding to reg¬ 
ular polygons (triangles, squares, pentagons, and 
hexagons) are used in the simulations apart from 
the circular grains. All grains have the same ma¬ 
terial density. Each simulation contains only one 
type of grain and these are monosized. Figure 
shows some snapshots of systems. The restitution 
coefficient is set in all cases to e = 0.01 and the 
friction coefficient n is varied between 0.1 and 2.0. 
The significantly low restitution coefficient allows 
for a rapid dissipation of the kinetic energy during 
the initial filling of the silo. The material density p 
of the grains is set to 0.01 Kg/m^ 

After the aperture is opened, we allow for any 
transient to fade before we calculate the flow rate 
and other properties. For any simulation, 15 inde¬ 
pendent realizations with different initial particle 
positions are carried out to estimate error bars. 

III. Results and discussion 

i. Mass flow rate 

We first consider the mass flow rate Q of grains as a 
function of the aperture size D for particles of same 
mass. Figure shows Q for different grain shapes. 
Since simulations are in 2D, the grains have been 
set with the same material density and area to yield 
the same mass m. Gircular grains present values of 
Q between 30% and 90% above those of polygonal 
grains. Squares and hexagons show larger fluctua¬ 
tion between independent realizations in compari¬ 


son with other polygons and disks. In particular, 
squares, which show the lowest flow rate, tend to 
clog rather often and data on Q for orifices < Sd 
could not be collected. Overall, the flow rate in¬ 
creases in the following order: squares, pentagons, 
triangles, hexagons and disks. 

We have fitted Q using the 2D Beverloo’s rule 

Q = Cp^{D-kdf/\ (1) 

Where p is the material density of the grains, and 
C and k are fitting parameters. Notice that the 
bulk density pt is generally used instead of the ma¬ 
terial density p which does not take into account 
the voidage of the system. In this case, we use p 
and let C to carry this information to avoid taking 
an arbitrary decision on the region where ph has to 
be measured. Table I shows the values of C and k 
obtained for each particle shape and Fig. [^includes 
the fitting curves. Disks, as well as polygons com¬ 
ply with the 3/2 exponent rule, at least within this 
narrow range of D explored. The values of C ob¬ 
tained confirm the general trend that the mass flow 
rate increases with shape in the order described 
in the previous paragraph. The fc-value indicates 
the intercept with zero flow rate and constitute an 
extrapolated estimate of the minimum size of the 
orifice necessary to observe flow. In practice, this 
estimate obtained with data at relatively large val¬ 
ues of D is little reliable. It has been shown that a 
more careful analysis of the flow at smaller D where 
intermittent flow can still be analyzed, suggest that 
a different equation to that of Beverloo is more ap¬ 
propriate by considering that only at fc = 1 would 
Q become strictly zero [1]. Nevertheless, the value 
of k gives an indication of how soon the flow rate 
drops as D is decreased. For disks, we find k « 2.0, 
which is consistent with experiments where a simi¬ 
lar range of D were used [S] (if smaller D are con¬ 
sidered k may drop down to « 1 [H [S]). Table I 
indicates that disks and pentagons are predicted to 
stop flowing at D « 2d. However, triangles and 
hexagons would only flow for somewhat larger D. 
Again, these are only rough indicators, since mea¬ 
surements at lower values of D may lead to different 
estimates, particularly if states of intermittent flow 
are included in the analysis. 

These results show that shape matters in granu¬ 
lar flow as shown by others [TTl [TH |T31 HU HU HS] • 
However, the variability of Q between different 
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Figure 2: Mass flow rate Q as a function of the 
aperture size D for different particle shapes with 
fi = 0.5. All particle shapes have the same area 
(mass). The continuous lines are fits to the 3/2- 
power oi D — kd (see fitting parameters in Table 

I)- 

faceted grains is rather limited in comparison with 
the significantly higher mass flow rate of disks. In 
general, setting apart triangles, the more vertices 
a polygon has, the higher Q is. This is consis¬ 
tent with the idea that polygons with many facets 
should tend to the behavior of circular grains. 

One may speculate that the use of faceted par¬ 
ticles changes particularly the frictional character 
of the interactions. While disks may roll on their 
contacts, polygons will mostly slide side-to-side on 
each other. A polygon may be forced to “roll” 
(rather than slide) on top of flat inclined surface 
if the friction angle 9 (i.e., 6 = tanp.) exceeds the 
maximum angle of stability Bstab (defined as the 
angle of inclination of a flat supporting base at 
which the vertical line that passes through the cen¬ 
ter of mass of the polygon meets one of its vertexes). 
Hence, it is expected that changing the friction co¬ 
efficient may lead to different effects on the poly¬ 
gons dynamics and hence on the flow rate. The 
corresponding values of beyond which a polygon 
would roll before sliding are > VS « 1.732 for 
triangles, /i. > 1 for squares, and > \/5 + « 

0.727 for pentagons. 

Figure shows Q for different fi for each grain 
shape. As we can see, friction does not induce any 
significant change in flow rate if fi > 1.0. How¬ 


Table 1: Fitting parameters for the Beverloo equa¬ 
tion in 2D as defined in Eq. Q. For squares, the 
poor statistics due to clogging events prevent us 
from extracting reliable fits. 


Shape 

C 

k 

Disks 

0.92 

2.01 

Triangles 

0.70 

2.30 

Squares 

- 

- 

Pentagons 

0.61 

2.09 

Hexagons 

0.82 

2.77 


ever, a considerable increase in Q is observed if we 
lower fj, below 0.5. The one exception is found in 
squares that are little affected by any change in /i. 
The increase in Q with decreasing /i is consistent 
with previous simulations of disks [8] . Interestingly, 
there are not clear evidences of a qualitative change 
when /i is higher than the maximum stability angle 
for each polygon. We presume that the interac¬ 
tion with other touching grains is as important as 
the interaction with any eventual grain that serves 
as primary support for a given polygon. Since the 
stability angle defined is only based on the sole in¬ 
teraction with a base, the expectation that a qual¬ 
itative feature may be observed at these /i values 
may be unjustified due to the complex multiparticle 
interactions. 

ii. Velocity, density and energy profiles 

The reduced flow rate of polygonal grains may be 
caused by a drop in the packing fraction (/> and/or 
a fall in the actual velocity of the grains. The 
packing fraction of polygonal grains in a rectan¬ 
gular box subjected to tapping has been previously 
considered [20]. It was shown that regular poly¬ 
gons able to tessellate the plane (triangles, squares 
and hexagons) in general arrange in higher pack¬ 
ing fractions as it may be expected. Pentagons, 
on the other hand, yield lower packing fractions 
in comparison with monosized disks. Despite this, 
the packing fraction during flow does not need to 
be connected with the values found in static sys¬ 
tems. Flow implies shear, which induces Reynolds 
dilatancy, which in turn reduces if the initial com- 
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Figure 4: Mean packing fraction distribution for 
disks (a) and triangles (b) averaged over a dis¬ 
charge. The aperture is D = lOd and the friction 
coefficient ^ = 0.5. Disks and triangles have the 
same area (mass). 


Figure 3: Mass flow rate as a function of D for 
different friction coefficients for each particle shape: 
disks, triangles, squares, pentagons and hexagons. 
All particle have the same area (mass) irrespective 
of their shape. 

paction is above the dilatancy onset. 

We have measured (j) profiles across the silo dur¬ 
ing a discharge for the case of disks and triangles 
of the same mass (see Fig. |^. The calculations 
are done by using the package Shapely for polygon 
intersections m- The time average of (j) during 
the discharge measured in cells of Id x Id are con¬ 
sidered. As we can see, 4> is lower for triangles 
than it is for disks in general. This is in contrast 
with the higher packing fractions found for trian¬ 
gles in tapped static systems [20]. It is clear that 
the reduction in (j) observed in the area of the outlet 
stretches farther into the silo for triangles. At the 
aperture, (j) is smaller overall for triangles and the 
effective “empty annulus” (i.e., the depleted region 
next to the edges of the aperture) looks thicker for 
triangles too. This explains the high values of k 
found in the fits of the Beverloo rule for triangles 
(see Table I), which is consistent with experimental 
results on sharp grains like some sand samples 
The reduced mass flow rate observed for polygons 
may be then caused by the reduction in (j). How¬ 
ever, the reduction with respect to circular parti¬ 
cles is less than 15%, which is insufficient to explain 
the at least 30% change in Q. Hence, one should 
expect that also the mean velocity of the exiting 
grains should be higher for disks. 

Figure shows the profiles for the vertical and 


horizontal velocity of disks and triangles. The mean 
horizontal velocity is different from zero only close 
the aperture in a V-shaped region consistent with 
convergent streamlines towards the outlet. Trian¬ 
gles show a more stretched region of non-zero hori¬ 
zontal velocity in agreement with the stretched re¬ 
gion of low (j) [see Fig. Bb)]. Disks present two 
hotspots above each edge of the outlet with high 
horizontal mobility in clear contrast with the tri¬ 
angles that present a smoother profile. From the 
mean vertical velocity one can observe clearly a 
more significant stagnant zone for triangles on each 
side of the aperture. This is consistent with results 
of several authors (see for example, (H] [14]) and 
with the intuition that sharp edged grains will dis¬ 
play a higher effective angle of repose. This taller 
stagnant zone seems to be responsible for inducing 
shear on a larger area for triangles, leading to the 
reduced (j). Also the vertical velocity at the aper¬ 
ture is in average higher for disks; roughly 9%. The 
combined effect of reduced 4> and reduced vertical 
velocity adds to a « 25% change in Q when com¬ 
paring disks and triangles at this opening width. 

In Fig. I^we show the translational kt and rota¬ 
tional kr kinetic energy profile for disks and trian¬ 
gles. kt is dominated by the vertical component of 
the velocity and therefore the features observed in 
Fig. E c) and Fig. [^d) are reproduced. An inter¬ 
esting feature observed in kr is that disks present a 
wide region over which grains rotate. This region 
extends all along the orifice width (and beyond, 
over the flat bottom base) at a height about 2d 
above the aperture. Triangles, in marked contrast. 
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Figure 5: Mean particle velocity: (a) Vx for disks, 
(b) Vx for triangles, (c) Vy for disks and (d) Vy for 
triangles. The aperture is D — lOd and the friction 
coefficient /i = 0.5. Disks and triangles have the 
same area (mass). The color scale indicates the 
velocity in ms“^. 

only acquire rotational energies less than half those 
of disks and over a limited region next to each edge 
of the aperture. This might be due to the frustra¬ 
tion of any initial spin a particle may acquire since 
collisions of its sharp edges with neighboring grains 
will tend to stop any rotational motion. 

At first sight, the reduced rotational energy of 
triangles seems counter-intuitive. If triangles have 
a reduced translational energy, one should expect 
that part of the energy input provided by the work 
done by the force of gravity might go into the ro¬ 
tational degrees of freedom. Notice however that 
the energy is also dissipated by collisions and fric¬ 
tion. It is clear then that triangles provide a much 
efficient dissipation of the energy for all degrees of 
freedom. 

iii. The effective size of a grain 

It is important to notice that comparing the be¬ 
havior of grains of different shapes conveys an arbi¬ 
trary decision on the parameter that decides when 
two grains are considered of the same size. In the 
previous section we have chosen to define as par- 


Figure 6: Mean kinetic energy of the particles: (a) 
translational for disks, (b) translational for trian¬ 
gles, (c) rotational for disks and (d) rotational for 
triangles. The aperture is D = lOd and the fric¬ 
tion coefficient /i = 0.5. Disks and triangles have 
the same area (mass). The color scale indicates the 
energy in J. 

tides of the same size those that have the same 
mass, irrespective of the shape. However, in some 
applications, one may be interested in comparing 
the flow of different polygons that have the same 
radius, side or perimeter. 

Figure]^ shows the mass flow rate for triangular 
particles of different sizes chosen to match the ra¬ 
dius, perimeter or area of the reference disks. We 
have also included data for triangles whose side is 
the same as the diameter of the disks. As we can 
see, since the actual radius of the triangle depends 
on the parameter used to define if the polygon is 
“of the same size as the disk”, the mass flow rate 
varies up to 25%. 

The importance of the arbitrary definition of the 
effective size of the grains when comparing different 
shapes was previously recognized El- However, 
the extent of the effect that this may cause in the 
mass flow rate was not studied. The 25% maximum 
difference observed in Q for triangles show that this 
cannot be neglected, although other shapes may be 
less affected. 

Notice that it is usual to compare particles of 
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Figure 7: Mass flow rate as a function of the 
aperture size D for disks and triangles with /r = 
0.5. Different triangles, that match the diameter, 
perimeter and area of the disks as well as having the 
side equal to the diameter of a disk, are considered. 


# # Disks (radius = d/2) 
i i Triangles (radius = d/2) 

^ Triangles (radius = 1.154 d/2, side = c(/2) 

^ ^ Triangles (radius = 1.21 d/2, perimeter = TTii) 
^ ^ Triangles (radius = 1.554 d/2, area = 7r(d/2)^) 
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different shape that have the same mass. However, 
there exist several works where authors use a char¬ 
acteristic length instead of the mass [in na E]. 
In some applications, due to design, economic cost, 
practical constraints, etc., fabricated objects of dif¬ 
ferent shape may or may not have the same mass. 
Estimations of the relative mass flow rate of such 
objects delivered by hoppers in comparison with 
circular counterparts would need to take into ac¬ 
count the relative size. 
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Figure 8: (a) Particle flow rate as a function of the 
aperture size D for different particle shapes with 
= 0.5 and the same particle radius, (b) Particle 
flow rate as a function of the aperture size D for 
different particle shapes with ^ = 0.5 and the same 
particle perimeter. 


iv. Particle flow rate 

It is worth mentioning that in some cases, the num¬ 
ber of grains flowed per unit time may be the im¬ 
portant parameter rather than the mass flow rate. 
Although this is proportional to the mass flow rate, 
the proportionality constant depends on the size 
of the grains and the mean </>■ Hence, depending 
on whether polygons are compared by their radius, 
side or perimeter, the actual particle flow rate, Qp 
varies. Naturally, if polygonal grains are compared 
by their area in 2D, the particle flow rate yields the 
same results as in Fig. [^simply scaled by the same 
factor (the particle density) for all shapes. How¬ 
ever, this is not the case if grains have, for example, 
the same diameter. 

In Fig. we present Qp for different polygons 


if the radius (a) or the perimeter (b) is used as 
the reference to compare sizes. As we can see, the 
relative position of the curves changes significantly 
depending on the parameter used to compare the 
particle sizes since the relative masses of the grains 
are different. If particles have the same radius [Fig. 
I^a)], pentagons and hexagons yield a slightly lower 
particle flow rate than disks, while squares deliver a 
larger number of particles per unit time. Triangles 
present up to 100% higher particle flow rates than 
disks. This is in dramatic contrast with the mass 
flow rate observed in Fig. If we use the perime¬ 
ter to compare grain sizes, most polygons have a 
lower particle flow rate than disks while triangles 
still display a larger flow in number of grains per 
seconds. 
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The results of Fig. |^show that the particle flow 
rate is particularly sensitive to the effective size of 
the grains. Squares, for example, may be delivered 
at a higher or lower rate than disks, in terms of 
number of particles, depending on whether they 
have the same radius or the same perimeter as 
disks. However, the mass flow rate of squares is 
lower than that of disks in both cases. 

IV. Conclusions 

We have considered the discharge of circular and 
polygonal grains from a flat bottomed silo. We 
have seen, in agreement with other workers, that 
for grains of the same mass, the mass flow rate 
is markedly lower for polygons in comparison with 
disks. We have found that this is caused by a com¬ 
bination of reduced mean packing fraction for the 
flowing polygons and also a reduced mean vertical 
velocity. The rotational kinetic energy of the poly¬ 
gons is also significantly smaller than that of disks. 

When comparing results for different shapes, it 
is necessary to arbitrarily select a parameter to de¬ 
fine when two grains can be said to be of equal 
size. We found that this selection, that may be 
constrained by specific applications, may lead to 
significant changes in the actual mass flow rates. 

We have discussed that the ratio of the parti¬ 
cle flow rate to the mass flow rate is a non-trivial 
function of the grain shape since it does not only 
depend on the mass of the grains but also on the 
actual packing fraction. As a consequence, it is dif¬ 
ficult to predict the particle flow rate from the sole 
knowledge of the mass flow rate. We have shown 
that the particle flow rate for a given shape, which 
can be the measure of interest in many industrial 
applications, can result higher or lower than for a 
different particle shape depending on the param¬ 
eter used to define their effective size. Hence, in 
a specihc application, in order to deliver the same 
number of grains per unit time for a different par¬ 
ticle shape, the aperture must be either reduced or 
widened depending on the parameter used to com¬ 
pare their sizes (e.g., diameter, mass, perimeter, 
etc). 
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